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ABSTRACT 

o : 
o . 

CN ' We describe the first results of a programme to obtain rotation curves of 2 ~ 1 

1-^ ■ disc galaxies in the near-infrared using the Ha emission line in order to study the 

pLn ' Tully-Fisher relationship. To put any observed evolution into perspective and to 

Q\ • investigate any possible selection biases, we constructed a control sample of low 

^ ■ redshift galaxies that had rotation velocities and images available for measuring 

^ ■ their dynamical, photometric, and morphological properties. Compared to local 

^ ■ objects with isophotal sizes similar to the high redshift targets, we find that 
our sample of galaxies with spatially resolved rotation curves, the most distant 

^ ■ sample so far {{z) ^ 0.9), clearly reveals a brightening of ~ 1.1 mag in the 

^ ' rest-frame B-band. The observed offset can be explained by a combination of 

O ■ increasing surface brightness, decreasing rotation speeds, and slightly smaller 

• disc scale lengths of the high redshift galaxies. 

^ ■ Subject headings: galaxies: evolution — galaxies: kinematics and dynamics — 

^ ■ galaxies: spiral 

> ' 

\-i ' 1. Introduction 

Like the fundamental plane for early type galaxies (Dressier et al. 1987) the Tully-Fisher 
relation for disc galaxies (Tully and Fisher 1977) embodies fundamental implications for 
the relationship between the mass of the galaxy, its star-formation history, specific angular 
momentum and dark matter content and distribution. 
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Broadly speaking, there are two competing models to explain the TuUy-Fisher relation. 
The first of these is that it is a consequence of self-regulated star formation in discs with 
different masses (e.g., Silk 1997), i.e., the competition of disc instability (which promotes star 
formation) with supernovae induced porosity (which inhibits star-formation). The model is 
not the complete answer, however, since it does not explain the mass-to-light ratios or 
the scale-lengths of the discs. In the second model the TuUy-Fisher relation is a direct 
consequence of the cosmological equivalence between mass and circular velocity (e.g.. Mo 
et al. 1998; Steinmetz and Navarro 1999). This formalism is part of what has become the 
standard model for the growth of structure - the hierarchical merging model in which the 
gravitational effects of dark matter drive the evolution of galaxies and large-scale structure 
(e.g., Kauffmann et al. 1993). Models of this type have the advantage of providing testable 
predictions about the sizes, surface densities, and rotation curves of galaxies as a function 
of redshift. However, as emphasized by Steinmetz and Navarro (1999), although the Tully- 
Fisher relation can naturally be explained by hierarchical merging models, the normalization 
and evolution of the TuUy-Fisher relation depend strongly on the prescription used for the 
star formation and on the cosmological parameters. 

It is now well established that massive disc galaxies exist out to redshifts z ~ 1 (Vogt 
et al. 1996, 1997; Schade et al. 1996; Lilly et al. 1998). For a given size scale, the number 
density of these disc galaxies is approximately the same at 2; ~ 0.8 as is observed locally. 
Overall, the results at moderate redshift {z ~ 0.5) are rather mixed. Depending on the 
sample selection, the technique used to estimate the rotation speed, the median redshift 
of the sample, and the wavelength at which comparisons are made, there are claims in 
the literature that the TuUy-Fisher relation either brightens or dims with redshift (see e.g., 
Vogt et al. 1996, 1997; Simard and Pritchet 1998; Rix et al. 1997; Ziegler et al. 2002). To 
help resolve this controversy and to push to higher redshift where evidence for evolution of 
disc galaxies will likely be stronger and more obvious, we set out to obtain spatially resolved 
rotation curves at the highest redshifts where large samples are available. Obtaining spatially 
resolved rotation curves becomes quite difficult at such high redshifts since [OII]A3727 is 
redshifted into a region with many strong night sky lines and the sensitivity of CCDs declines 
rapidly. Instead we have chosen to observe the Ha line, reshifted to the infrared J-band. Ra 
is typically 2.5 times stronger than [OH] (Kennicutt 1992) and being at longer wavelength, 
is less affected by dust extinction. 
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2. Sample Selection 

To gauge the evolution of the Tully-Fisher relation we compare a local with a high 
redshift sample of highly inclined, "normal" spiral galaxies. The high redshift sample con- 
sists of objects with measured spectroscopic redshifts. To be able to measure Ha in the 
near-infrared we required z > 0.6. Targets were selected from the CFRS / LDSS redshift 
surveys, the clusters MS1054, AC103 and the Hawaii deep field SSA22 (Lilly et al. 1995; 
Glazebrook et al. 1995; van Dokkum 1999; Couch et al. 1998; Cowie et al. 1996, respectively, 
and references therein). Furthermore, we included two targets from the VLT science archive. 
For the majority of these sources we have HST I-band^ or at least VLT J-band images. For 
all of them we have obtained Ha ISAAC long-slit spectra to determine spatially resolved 
rotation curves. For this purpose we have selected these sources to be spatially extended on 
the images (diameters mainly >2") and to have high inclination angles {i > 30°, (^) ~ 60°). 
Furthermore, we chose only objects with exponential light profiles to ensure that we were 
observing disc galaxies. The mean redshift of our high-z sample is z ^ 0.92 (see Tab. 1). 

To construct a complementary sample of local galaxies we searched the literature for ob- 
jects with accessible CCD B-band images (obtained from the NASA Extragalactic Database, 
NED, or ESO science archive), pubhshed distances based on primary distance indicators or 
mean group velocities (Freedman et al. 2001; Pierce and TuUy 1992; Puche and Carignan 
1988; TuUy and Pierce 2000; Ferrarese et al. 2000; Giuricin et al. 2000; Drozdovsky and 
Karachentsev 2000; Karachentsev and Sharina 1997; Karachentsev et al. 2000; Distefano 
et al. 1990; Abell et al. 1989; Duval and Monnet 1985, LEDA database^) and an estimate 
for the rotation speed based on either an HI or a resolved Ha rotation curve (Prugniel et al. 
1998, and references therein). Total magnitudes for all local sample galaxies where ob- 
tained from RC3 (de Vaucouleurs et al. 1991) or LEDA. As a consistency check, the object 
magnitudes, sizes, rotation speeds and the resulting Tully-Fisher relation were compared 
to the RC3 catalogue and the data from TuUy and Pierce (2000). From this comparison, 
we conclude that our sample of low redshift galaxies is in good agreement with a random 
subsample of a magnitude limited complete sample from the RC3. 



^Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data 
archive at the Space Telescope Institute. STScI is operated by the association of Universities for Research 
in Astronomy, Inc. under the NASA contract NAS 5-26555. 

^The Lyon Meudon extragalactic database (http://leda.univ-lyonl.fr) 
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3. Observations 

We obtained J- and H-band spectroscopy during four observing runs with the faciUty 
near-infrared spectrometer ISAAC at the VLT (Moorwood 1997; Moorwood et al. 1998). 
We used the medium resolution spectroscopic mode with a slit-width resolution product of 
R ^ 3000 - 5000 and a spatial pixel scale of 0.1484" . 

In ESO-period 65 and 66 (April - September 2000 and October 2000 - March 2001 
respectively) we obtained the first set of data in service mode. With total integration times 
of 3 hours and a 0.6" -slit {R 5200), which we rotated according to the position angle of 
the objects, we could successfully observe four targets. The target selection and successive 
estimation of position angles for this set we performed on 1" resolution images acquired in 
the J-band with ISAAC (20min integration time) before the spectroscopic observations. In 
addition to these J-band images, we obtained I-band images taken with the Canada Prance 
Hawaii Telescope (CFHT) from that telescope's archive. 

Furthermore, we conducted three other runs in visitor mode in February and September 
2001 and September 2002. This time the set-up included a l"-sht {R « 3100) and integration 
times varied, depending on the Ha flux of the targets, between 1 and 3 hours. Moreover, we 
selected targets on the basis of HST/WFPC2 F814W images, which we also used to orient 
the slit along the galaxy major axis. Since most of the targets were too faint to be acquired 
directly we used a blind off-set from a "bright" star, calculated from the HST images. Using 
this prescription we observed another 17 objects successfully. The seeing during these four 
runs ranged from 0.4" to 0.8". 

4. Data Reduction and Analysis 

To reduce the ISAAC spectra we used standard data reduction packages (ESQ eclipse & 
IRAF). The data were flat fielded and dark-subtracted. After dark-subtraction and cosmic 
ray rejection we checked all frames individually for shifts in the wavelength direction due to 
flexure and applied a reverse shift if necessary. Before combining the frames we calculated 
a wavelength solution and distortion correction. Finally, we subtracted a background and 
smoothed the spectra with a Gaussian of FWHM approximately equal to the number of 
pixels per resolution element. 

Firstly, we extracted one dimensional Ha-profiles across the spectral and spatial regions 
of the array that contained Ha-flux. Typically, even with the increased signal-to- noise gained 
in the continuum by summing up the flux in this way, we could not detect the continuum 
significantly. We then fitted the one dimensional Ha-profiles (typically assuming no galaxy 
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continuum) by one or two (in the case of a double horned profile) Gaussians. At the level 
where 20% of the total integrated flux was reached we measured the width of the profile. 
Applying the corrections given by Tully and Fouque (1985) we obtained a first estimate for 
the terminal rotation velocity of our sample galaxies. 

In all cases but one, however, we could derive a terminal rotation velocity from a resolved 
rotation curve model. To extract the data points for the model rotation curve we developed 
a special adaptive curve tracing algorithm. We then fitted the resulting position-velocity 
diagram with a simple model of a fiat rotation curve (which is appropriate for not too-late 
disc galaxies). This step-function we convolved with a Gaussian with a full- width at half 
maximum of the seeing. Moreover, the convolution included weighting with an exponential 
function with a scale-length as measured from our disc surface brightness fits on the HST 
or VLT images assuming that line and continuum have the same light distribution. Tests 
show that our modeling delivers robust rotation velocities, which are not very sensitive to 
the exact weighting scheme, the impact of nuclear bulges, the signal-to-noise ratio and in 
particular, how well the observationally "flat part" of the rotation curve is sampled by the 
data. We have measured the seeing from the J-band acquisition images taken before the 
observations, the acquisition images of a subsequent target, and the values from the visible 
seeing monitor. Given our long integration times and the limited sampling of the seeing at 
the observing wavelength, the seeing is somewhat uncertain. We have therefore let it be a free 
parameter in the fitting process. The values obtained from an unconstrained fit correlate well 
with the seeing values estimated from the acquisition images and seeing monitor, though. 
Some of the final fits to the data are shown in Fig. 1. Although only 50% — 60% of the data 
show a clear turn-over or extend out to 2.2 disc scale lengths (which is were the rotation 
curve of a self-gravitating thin disc reaches its maximum, see Binney and Tremaine 1987) 
the remaining data do not exclusively consist of slow rotators only. In fact, the distributions 
of rotation speeds of high and low quality data can not be distinguished from each other (see 
Barden et al. 2003). Furthermore, the rotation speeds derived from the two dimensional 
data are not significantly smaller than the ones from the one dimensional spectra, which 
have much higher signal-to- noise (see Barden et al. 2003). Moreover, according to the 
parametric fits to over 2000 local disc galaxies of GiovaneUi and Haynes (2002) at radii 
varying from 1...1.5 x one should already measure at least ~ 80% of the asymptotic 
rotation velocity. Since the majority of our high galaxies extend further out we estimate 
rotation velocities within < 20%. Therefore, we are confident that we do not underestimate 
the rotation velocities significantly. Finally, the rotation velocities have to be corrected for 
the inclination of the galaxies. We have derived inclinations from surface brightness profile 
fits to our HST or ISAAC images. 

To estimate rest-frame B-band magnitudes we compiled all available photometric in- 
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formation for our sources (two to five independent measurements). For each magnitude we 
calculated the corresponding rest-frame wavelength and corrected the magnitudes for galac- 
tic foreground extinction (Schlegel et al. 1998). Then we fitted the extinction corrected 
magnitudes and corresponding rest-wavelengths with spectra from the atlas of Kinney et al. 
(1996). We note that our final rest-frame B-band magnitudes are relatively insensitive to 
the models used since all of the galaxies had a photometric band available that was very 
near to the B-band rest wavelength. Finally, we corrected the resulting rest-frame B-band 
magnitudes for inclination dependent internal extinction according to the method of TuUy 
and Fouque (1985) and determined a face-on B-magnitude (Table 1). 

The parameters for the local comparison sample were determined and corrected in ex- 
actly the same manner as for the high redshift galaxies. Rotation speeds were extracted 
from HI or optical rotation curves. Only one object was included where a spatially-resolved 
rotation curve was not available. Comparison of the resolved data with the HI line width 
from RC3 usually agree within ±20 km s~^ (see Harden et al. 2003). We calculated inclina- 
tion corrections based on inclinations from our own profile fits in agreement with the values 
from RC3 or LEDA. Total magnitudes were obtained directly from RC3 or LEDA and were 
corrected for galactic foreground extinction and internal extinction. Since differences in the 
isophotal extent between the images of the local and high redshift galaxies can be quite 
dramatic (up to two magnitudes), the light profiles of the local galaxy sample were fit only 
out to a characteristic limiting radius, which was approximately the average isophotal depth 
of the images of the high redshift galaxies (taking into account cosmic surface brightness 
dimming; see Barden et al. 2003). 



5. The Tully-Fisher Relation at z ^ 0.9 

Comparing the isophotal-size-cuUed local (our comparison sample) and high redshift 
galaxy samples (see Barden et al. 2003) we find an offset for the high-z data from the 
Tully-Fisher relation of AM = —1.1 ± 0.2 magnitudes in a A-dominated universe {Hq — 70, 
JIm = 0.3, Qa — 0.7; Fig. 2). This "brightening" of the distant galaxies is significant 

and holds for any reasonable cosmology (AM = —0.5 ± 0.2 for VIm = 1 and Qa = 0; 
AM = —1.0 ± 0.2 for Qm = 0.05 and Qa = 0). The limited range in isophotal size makes 
measuring a slope of the Tully-Fisher relation impossible. 
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6. Discussion and Conclusions 

The observed offset in the Tully-Fisher relation at {z) ^ 0.9 for large galaxies is partially 
a consequence of the fact that we observe an increase of central surface brightness of 1.7 mag 
arcsec"^ and a slight decrease of ~ 30% in the disc scale lengths of the high redshift sample 
compared to the local control sample (see Barden et al. 2003). Using the relationship 
between disc magnitude, central surface brightness, and disc scale length, namely, oc 
—5 log Rd — 2.5 log {1 + z) + fj,o — const, we find that the differences in average central surface 
brightness and average disc scale lengths imply an offset in the zero-point of the T-F relation 
between the two samples of about —0.3 ± 0.4 magnitudes. In addition, by comparing the 
rotation speeds of the high redshift galaxies and the local comparison sample, we find an 
average offset of Avrot = 50 ± 15 km s^^ of the high redshift sample (see Barden et al. 
2003). Thus, between z = 0.9 and the mean mass of galaxies has increased by a factor 
of 2.6 ± 0.4 (assuming mass oc f^^j, and estimating the total mass out to 2.2 disc scale 
lengths measured from fitting the light profile). Therefore, if we use a parameterization of 
the Tully-Fisher relation as given by TuUy and Pierce (2000), adjusting the zero-points by 
-0.3 mag to account for the average differences in surface brightness and disc scale length 
and decreasing the average rotation speed by 25% yields an offset from the local Tully-Fisher 
relation consistent with that observed (AM — — 1.1 ± 0.5 mag). 

Having the highest redshift sample puts us in a unique position to discuss the evolution 
of the Tully-Fisher relation from a new perspective. All of the studies to date have found a 
significantly smaller amount of evolution compared to the results presented here. The most 
obvious way to reconcile this difference is that the offset in the Tully-Fisher relation grows 
with increasing redshift. The studies with the highest median redshift currently are Ziegler 
et al. (2002) and Vogt et al. (1996) with a mean redshift of (z) ~ 0.4 — 0.5, compared with 
our mean redshift of {z) fti 0.9. The time difference is therefore about 2-3 Gyrs for any 
reasonable cosmology, clearly sufficient for galaxies to undergo significant evolution. 

Unfortunately, comparing these results with model predictions is not very straight- 
forward since the exact offset predicted depends on the assumed physics and phenomenology, 
such as the details of star- formation (e.g., Steinmetz and Navarro 1999; Mo et al. 1998). 
Recently, within the frame work of galaxy evolution (e.g.. White and Rees 1978), a number 
of authors have made predictions about the general evolution of the galaxy population. 
For example, Nagamine et al. (2001) predict that from 2; = to 1, using the cosmological 
parameters we have assumed here, that the overall luminosity function of galaxies in the 
B-band will brighten by about 1.15 mag and that the characteristic mass of galaxies will 
decline by about a factor of 2.4. Assuming that the galaxies selected for study here are 
a subset of this population, our numbers are in reasonable agreement with the predicted 
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values. Unfortunately, the resolution of the simulations by Nagamine et al. (2001) are too 
low to check whether the sizes of discs evolve self-consistently in this context. 

We conclude that large disc galaxies have undergone substantial evolution from 2; = 0.9 
to the present. The offset in the TuUy-Fisher relation is about —1.1 mag in the rest- frame B- 
band, and is related to an overall increase in the central surface brightnesses of disc galaxies 
and a decrease in rotation speed. We also find that the slope of the TuUy-Fisher relation, due 
to size constraints that inevitably result from needing to obtain spatially-resolved rotation 
curves, cannot meaningfully be estimated. In a subsequent paper (Bardcn et al. 2003), we 
will discuss in detail the change in rotation speeds, disc masses, disc scale lengths, and disc 
angular momenta of this sample of high redshift galaxies relative to the local comparison 
sample, and make quantitative comparisons with theoretical models of galaxy evolution. 

We wish to thank the ESO OPC for the generous allocations of telescope time and the 
Paranal staff for their expert support and assistance during the observations. 
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Table 1. The high redshift sample 
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^galaxy observed by White et al. (ESO proposal ID: 63.O-0372(A)) 

Note. — (1) galaxy identification (Lilly et al. 1995; Glazebrook et al. 1995; van Dokkum 1999; Couch 

et al. 1998; Cowic ct al. 1996) (2) redshift, (3) isophotal radius at a surface brightness of 25 mag arcscc"^, 
(4) ratio of major over minor half axis, (5) inclination, (6) internal extinction (TuUy and Fouquc 1985), 
(7) galactic extinction (Schlegel et al. 1998), (8) total rest-frame B-band magnitude, (9) total absolute 
rest-frame magnitude plus error (assuming Hq = 70., Q,m = 0.3, Ha = 0.7), (10) rotation speed corrected 
for inclination plus error 
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Fig. 1. — Typical 2-D spectra for some galaxies of our high redshift sample. The grey-scale 
shows the 2-D Ha spectra, the squares represent the adopted position-velocity measurements 
with the best-fit model rotation curve over-layed as a dashed line, and the vertical dashed 
and/or sohd lines show the isophotal radius R25 at a surface brightness of 25 mag arcsec"^ 
and 2.2 times the disc scale length R2.2 — 2.2 x i?^, respectively. 
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Fig. 2. — TuUy-Fisher relation for the high and low redshift samples. Circles correspond to 
high redshift galaxies, diamonds to local galaxies. Pluses mark isophotal sizes larger than the 
required minimum R25 > 12.5 kpc. Only local galaxies shown as large symbols are included 
in the control sample in order to match the RC3 statistics. The thin lines mark the local 
relation by TuUy and Pierce (2000) (thin dashed hne) and its iScr error bars (thin dotted 
hnes). A fit to the high redshift data results in a TuUy-Fisher relation -l.lmag offset from 
the local TuUy-Fisher (thick sohd line). Our reconstruction is plotted as the thick dashed 
line. The agreement between fit and our high redshift calibration is astonishing. 



